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Abstract 

The method of solution of the initial value constraints for pure 
canonical gravity in terms of Ashtekar's new canonical variables due 
to CD J (see (TJj) is further developed in the present paper. There are 
2 new main results : 

1) We extend the method of CD J to arbitrary matter-coupling 
again for non-degenerate metrics : the new feature is that the 'CDJ- 
matrix' adopts a nontrivial antisymmetric part when solving the vec- 
tor constraint and that the Klein- Gordon- field is used, instead of the 
symmetric part of the CDJ-matrix, in order to satisfy the scalar con- 
straint. 

2) The 2nd result is that one can solve the general initial value 
constraints for arbitrary matter coupling by a method which is com- 
pletely independent of that of CD J. It is shown how the Yang-Mills 
and gravitational Gauss constraints can be solved explicitely for the 
corresponding electric fields. The rest of the constraints can then be 
satisfied by using either scalar or spinor field momenta. 

This new trick might be of interest also for Yang-Mills theories on 
curved backgrounds. 
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1 Introduction 



The general solution of the initial value (iv) constraints with respect to the 
vector and the scalar constraint within the Ashtekar- framework (see f3j) for 
pure gravity has a strikingly simple structure when restricting it to non- 
degenerate metrics. 

The gravitational action in the new canonical variables introduced by Ashtekar 
(see ref. ||) has the form 

E S= f dtf d 3 x{P?< - [~4{ E Qi) + N a ( E V a ) - kN( e C)]} . (1.1) 

Here Pf = — i/2/ty / det(g)e" where k/(8tt) is Newton's coupling constant and 
e" is the triad of the metric q := (qw,) on the hypersurface S while u) l a is 
the Ashtekar-connection. We use indices a,b,c,.. from the beginning of the 
alphabet to denote valences of tensors defined on the initial data hypersurface 
£ while indices i,j,k,.. from the middle of the alphabet describe the 0(3)- 
gauge group structure of a generalized tensor. CD J (see |J) have shown that 
the ansatz 

P- = %jB] , (1.2) 

where Bf := \e abc VL\ c is the magnetic field with respect to the Ashtekar- 
connection, Q l ab being its field strength and e abc is the totally skew (metric- 
independent) tensor density of weight 1, solves both the gravitational part 
of the vector constraint 

E V a := n*jj? = (1.3) 
and the gravitational part of scalar constraint 

E C := K b P«P b k e ijk = (1.4) 

provided the CDJ-matrix \l/ is subject to the following conditions : 

* T = * and (tr(^)) 2 - tr(* 2 ) = . (1.5) 

Here A T means the transpose of the matrix A. Inserting this ansatz into 
the gravitational part of the Gauss-constraint one obtains when using the 
Bianchi-identity (T> is the gauge covariant differential acting on arbitrary 
generalized tensors) 

E Qi '■= T^aPi — Bj'Da^ij = (1.6) 
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which is now a differential condition on \1/ and cannot be solved by purely 
algebraic methods any more. Up to now there are no solutions known to 
equation (1.6) for full gravity. 

Furthermore, it should be stressed that even if one had also the general 
solution to equation (1.6) then one would only have obtained the constraint 
surface of the 'non-degenerate sector' of the gravitational phase space, not 
its reduced phase space since one did not factor by the gauge orbits yet. 
Nevertheless, this method could be of some importance for obtaining the 
reduced phase space for full gravity as might be indicated by the fact that it 
is of some help in model systems (see j4j). 

In this paper we are going to discuss the iv constraints in terms of Ashtekar's 
new variables for arbitrary matter coupling. The analysis will be in the 
canonical framework as in equation (1.1), so all the fields will be subject to 
4 different types of constraints, namely 

1) the gravitational Gauss constraint (Lagrange multiplier : —col) 

2) the Yang-Mills Gauss constraint (Lagrange multiplier : — A\) 

3) the vector constraint (Lagrange multiplier : N a , the shift vector) and 



4) the scalar constraint (Lagrange multiplier : N, where N := Jdet(q)N is 
the lapse function). 

The matter sector which we are going to discuss consists of (we will not dwell 
on how to derive the 3+1 form of the various actions, for details see ref. 0; 
the 3+1 form of the Higgs-action is not derived there but it can be obtained 
by a calculation similar to that for the Klein- Gordon action so we can omit 
this here) 

matter g = KG g + W g + C g + YM g + H g ? ^ 

i.e. a collection of (real) Klein- Gordon- fields with arbitrary potential V{4>) 
(we do not display the summation over the different scalar fields) 

KG S = I dtf rf 3 x{vr0 - [iV a 7T0 ia + ]-N{n 2 + det(g)(g a V,a0, fe + V (</>)))]} , 

(1-8) 

a (collection of) (complex valued rather than Grassmann- valued) Weyl spinor 
field(s) which couple only to the self-dual part of the spin-connection (this is 
explained in refs. |J and || ; we write down only one spinor field which may 
stand for an arbitrary number of Weyl-fields of possibly both chiralities), 
including an arbitrary spinor potential (e.g. the usual mass term when at 
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least 2 spinor fields of both chiralities are present) 
W S = f dt f d 3 x{7r T iP - [-u; t V T 7# 

J R JY1 



+N a n T V a ij - ±N{KPy T r{D a i\> + det(q)V(ij A , ir A ))}} , (1.9) 



where Tj = —i/2ai, <7i the usual Pauli-matrices, V a ip := [d a + oo l a Ti\ip and 
7r = (77,4), ^ = C0 A ), A =1,2 SU(2) spinor indices, 

a Yang-Mills-field for a semi-simple gauge group G (such that the Cartan- 
Killing metric (d/j), I,J,-- internal indices of this gauge group, is non- 
degenerate; if Tj are the generators of the Lie algebra LG of G then du = 
triTjTj), [Tj,Tj] = f u K T K ) 



YM S = I dt I d^xiE^Ai - [-A*V a E? 

JR JT. 



+N a F I ab E b I + £ Nq ab d IJ (E«E b j + flffl*)]} , (1.10) 

where Ef is the YM electric field, Bj := l/{2g 2 )d IJ e abc F b J c , F T ab being the YM- 
field strength, is the YM magnetic field and g the YM coupling constant, 
a cosmological constant (A) term 

C S= I dt I d 3 xNdet(q)A (1.11) 

J R -J S 

and a Higgs- field for either the gravitational gauge group 0(3) or the YM 
gauge group G 



L ' H S = f dt f d 3 x{iri<j) 1 - [-wieij k <jy>Tx k + N a n{D a ^ 

JR JT 
1 



YMH 



+ ^N^Wj + det( g )( g %(P a <^(W) + W)))]} or(1.12) 
S = f dt I d 3 x{7i I <p I - [-Mfu K (j)\ K + N a ^jV a ^ 

JR JT, 

+ ^N(d IJ m7Tj + det(g)(g a % J (P a / )(P 6 J ) + \/(/)))]}(1.13) 

After this brief introduction of our notation, we can outline the plan of the 
paper : 
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In section 2 we will show how the method of CDJ can be extended to 
include arbitrary physically relevant matter couplings. The deviations from 
the source-free case are that the CDJ-matrix \1/ fails to be symmetric and that 
the scalar constraint is not solved for \l/ but for one Klein-Gordon-momentum 
it. The equation for tt is purely algebraic and of 2nd or 4th order repectively 
depending on whether the YM field is coupled or not. This is just sufficient 
in order that the scalar constraint be solvable algebraically, for example, by 
the methods of Cardano and Ferrari (see 0). As for the source-free case, 
this method does not solve the gravitational Gauss constraint. 

In section 3 we show that by employing a method which is completely 
independent of the CDJ-framework it is possible to solve all constraints of 
general relativity coupled to arbitrary matter. An advantage of this method 
is that the solutions of the constraint equations are remarkably simpler than 
those obtained by the method presented in section 2. 

First we show how the gravitational Higgs-field can be used to solve both the 
gravitational Gauss constraint and the scalar constraint. The equations for 
the Higgs field are at most quadratic. The YM-Higgs-field can (in general) 
be used in order to satisfy the YM Gauss constraint. Now it is the vector 
constraint which is not solvable by purely algebraic methods for the gravita- 
tional field. However, by using again the gravitational Higgs field and at least 
3 Klein-Gordon fields one can solve both the vector and the scalar constraint. 
By a similar line of approach one can also solve only the YM Gauss-constraints 
purely for the YM electric field. The rest of the constraints can then be sat- 
isfied by purely algebraic methods if at least 4 Weyl spinor fields are present 
which, of course, corresponds to the physical reality. This latter option has 
the advantage that one can forget about scalar fields altogether which have 
not yet be proved to exist at all. 

The results of section 3 may be of interest also for the canonical approach to 
(pure) YM theory on curved or flat background metrics. 

The paper concludes with an appendix in which some explicit formulas 
for special cases of couplings, i.e. when the YM-field is absent, are given, 
valid when applying the method of section 2. There we also display the 
formulas of Cardano and Ferrari for the reader who wants to arrive at an 
explicit solution of the initial value constraints including a YM-field when 
applying the framework of section 2. 
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2 Solving the vector and the scalar constraint 

The full vector constraint reads 

V a = W ab P* + 7r0, a + rr T V a tlj + Fi h E\ + %{D a ^ + iv^ 1 . (2.1) 
As CD J (ref. [[I]]) we make the ansatz 

P- = %jB% (2.2) 

and obtain 

where B^B l b = 5 b , B?B J a = 5j. This ansatz poses no restrictions on Pf as 
long as the magnetic fields are non-degenerate (which however excludes flat 
space field configurations). 

The important step is to decompose the CDJ-matix into its symmetric and 
antisymmetric part 

^ := S + A, S T = S, A T = -A (2.3) 

because only A enters the vector constraint which then can in fact be solved 
for A : let 2£ k := e ijk A i:j , B := det(Bf), then 

& = ^Tnl^J + [H^ 7 ^ + F abEl + 7T t V a ^ + TTjVa^ 1 } =: 7T % + 6 k (2.4) 

Zt> Zjd 

where the quantities rjk, Ok do not depend on tc. Contracting equation (2.4) 
with e.ijk we obtain 

A : = 7rT + R where T T = —T, R T = —R , (2.5) 

i.e. the matrix A is linear in tt, homogenous only if the matter different from 
the KG- field is absent. 

We now insert this into the scalar constraint whose complete expression is 
given by 

C = -^P/P^ + If^+det^^a^ + V^))] 

-A{KP^ T n V a ^ + det(q)V(^ A , tta)) + A det(g) + ^-q ab d IJ (E«E b j + B^B)) 
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+ l -[5 ij ir i <K j + det(q)(5 t3 q ab (V a f)(V b ^) + V (ft))] 
+ ^[d IJ 7Cj7rj + det(g)(rf /J g afe (P a / )(W 7 ) + W))] 

= ■ E C+ KG C+ W C+ C C+ YM C+ EH C+ YMH C 

When inserting for Pf the various terms involve different powers of A which 
we now discuss in sequence : 

~- E C = e abc B<rB a m B b n e^^ kn = Be mni e jki ^ jm ^ kn = £?((tr(M/)) 2 -tr(M/ 2 )) 

(2.7) 

Since tr(A) = tr{A T ) = -tr(A) = 0, tr(SA) = tr(A T S T ) = -tr(AS) = 
—tr(SA) = for any symmetric matrix S and any antisymmetric matrix A, 
it follows easily that 

E C = -KB((tr(S)) 2 - tr(S 2 ) - tr(A 2 )) , (2.8) 

such that 7r enters E C only purely quadratically without a linear term. 
Let (fr^B? := v i: then 

det(g)g a V,a0, fe = -4 K 2 i^V,a0, fe , (2.9) 

but Pj a 5 a = (Sij + Rij)vj because T^Vj = eij k VjV k /(2B)(j) ja = 0. Hence 

det{q)q ab <f) a <f) b = -Ak 2 v t ((S+R) t (S+R))v = -4K 2 tr{{S 2 -R 2 +[S, R])v®v) 

(2.10) 

which is independent of ir. 
Now 

det(det(g)<f b ) = (det(g)) 3 (det(g))- 1 = (det(g)) 2 
= -(2 K ) 2 det(^) = -(2 K det(^)) 2 (2.11) 

from which follows (up to a sign) that 

det(g) = 2iKBdet{V) . (2.12) 

Since a term cubic in A enters det(^>) only through det(A) (this follows from 
the fact that \& = A+S and that the determinant is a totally skew multilinear 
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functional) which vanishes in 3 dimensions, det(q) is also only quadratic in 
ii. The explicit formula 1 is given by 

det(tf) = ^e ijk e lmn (A + S) u (A + S) jm (A + S) kn 

= ^\ e ^ keimn ^ il ^^ m ^ kn SSuSjmAkn + 3SuAj m A kn + AuAj m A kn ) 

= det(S) + det(A) + ^3\5 l [t 5f5 n k] (S tl S jm A kn + S u A jm A kn ) 

= det(S) + tr(SA 2 ) - \tr(S)tr(A 2 ) . (2.13) 

Altogether we have therefore 

KG C = \[k 2 - AK 2 tri(S 2 - R 2 + [S,R\)v®v) 

+2iKB(det{S) + tr(SA 2 ) - \tr(S)tr(A 2 ))V {<j)))} . (2.14) 

The Higgs-sector differs from the KG-sector as far as the appearence of ix is 
concerned only in that T^B^cf)^ = 0, while TijBjV a <f) k , TijBjV a <f) k ^ 0. 
Let v) := B*V a (f)\ v] := B^V^ 1 then 

"C = l -P>«* i -**S,M { S>-A' + lS,A\y 9 *) 

+2iKB(det(S)+triSA 2 ) - ^r(S)*r(A 2 ))y(^))] 
YMH C = \[d IJ Tx^j-AK 2 djjtr({S 2 - A 2 + ®v J ) 

+2mB(det{S)+tr(SA 2 ) - \tr{S)tr{A 2 ))V {(f) 1 ))] • (2-15) 
There is only a quadratic appearence of n again. 

The contribution by the YM-sector however gives rise to a quartic constraint 
for n : 
We have 

q ab = det(g)£X = j^e^e bef e^e imn P^P^pf 



1 This formula can also be obtained by applying the theorem of Hamilton-Cay ley as 
was pointed out to the author by Ted Jacobson in private communication 
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(2k) 4 B 2 

— 2 det(g) trst€uvw ^J^k ^ j^ks^mu^nvB a Bf 
= |^3!^^(* r *) rB (* T *) WJ B* a S 6 " 

= WTT^i( tr (^)) 2 ~ tr((* T *) 2 ) + 2(* r *) 2 _ 2tr(^)^] tw ^ a Bf) 
2det(g) 

Since tf T tf = (S — A)(S + A) = S 2 - A 2 + [S, A], tr([S, A]) = 0, we conclude 
that 

q ab = ^^[(tr(5 2 -A 2 )) 2 -tr((5 2 -A 2 + [5,A]) 2 ) + 2(5 2 -A 2 + [5,A]) 2 

-2tr(S 2 - A 2 )(S 2 - A 2 + [S^A^BlB™ . (2.17) 

As long as det(g) 7^ one is allowed to multiply the scalar constraint by 
det(g) so that it becomes 4th order in n because the contributions of the 
other fields are at most quadratic in n and det(\l/) is, according to formula 
(2.13), also only quadratic in ir. 

The final step is to collect all terms and to determine the coefficients of the 
various powers of n. We have 

det(tt) = [det(S) + tr(SR 2 ) - -tr(S)tr(R 2 )] + ir[tr(S(RT + TR)) 

2 

~tr(S)tr(TR+ RT)] + n 2 [tr(ST 2 ) - \tr(S)tr(T 2 )} 
=: (an 2 + bn + c)/(2inB) (2.18) 

and 

= {S 2 - R 2 - [S,R]} -n{RT + TR- [S,T]} - tt 2 T 2 , 
(^) 2 = [(S 2 -R 2 -[S,R]) 2 ] 

-7i[(RT + TR- [S, T])(S 2 -R 2 - [S, R}) + (S 2 - R 2 - [S, R}) 
(RT + TR- [S, T})] - tt 2 [(RT + TR- [S, T})T 2 
+T 2 (RT + TR- [S, T}) - (S 2 - R 2 - [S, R]) 2 ] 
+n 3 [(RT + TR- [S, T])T 2 + T 2 (RT + TR- [S, T])] + tt 4 T 4 
=: A 4 7i 4 + A 3 rr 3 + A 2 n 2 + A l7 v + A (2.19) 
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such that 

2det(g) b 
9a6 (2«) 4 £ 2 1 j 

= [(tr(S 2 - R 2 )) 2 - Antr{S 2 - R 2 )tr(RT) - 2Tr 2 (tr(S 2 - R 2 )tr(T 2 ) 

-2(tr(RT)) 2 ) + 4ir 3 tr(RT)tr(T 2 ) + ir 4 (tr(T 2 )) 2 - (tr(A 4 )ir 4 + tr(A 3 )ir 3 
+tr(A 2 )n 2 + tr(Ai)ir + tr(A )) + 2(A 4 vr 4 + A 3 tt 3 + A 2 rr 2 + A x n + A ) 
-2{tr{S 2 - R 2 ){S 2 -R 2 - [S, R}) - n(tr(S 2 - R 2 ){RT + TR- [S,T]) 
+2tr(RT)(S 2 — R 2 — [S, R])) 

-7t 2 (tr(S 2 - R 2 )T 2 + tr(T 2 )(S 2 - R 2 - [S, R]) - 2tr(RT)(S 2 - R 2 - [S, R])) 
W{2tr(RT)T 2 + tr(T 2 )(S 2 - R 2 - [S, R])) + 7rHr(T 2 )T 2 )] tj . (2.20) 

We first collect the powers of n contained in the non-YM part of the scalar 
constraint and then multiply by det(g). 

E C + KG C+ W C+ C C + EH c + YMH Q 

= {-KB((tr(S)) 2 - tr(S 2 ) - tr(R 2 ) - 2kHt{{S 2 - R 2 + [S, R])v ® v) 
+c^V(cf>) - 4(kB](S + R^nV^ + cV(^ A , tt a )) + cA 

+cV((f> 1 ))] + \[d IJ Tim.j - AK 2 d u tr([S 2 -R 2 + [S, RWv 1 ® v J ) + cVty 1 ))]} 

+ {2KBtr(RT) + h -V{<t>) - 4(kB?7> 3 \{D a ^ + bV(^ A , n A )) 
Z J 

+bA + hAK%tr([KT + TR- [S,T]]]v l ® v j ) + bV {<¥))} 
+ ^[AK 2 d L jtr{{RT + TR- [S, T^v 1 ® v J ) + bV(^))]}n 
+ {nBtr{T 2 ) + + aV{4>)\ - AaV(^ A , ^ A ) + aA 
+ -[4K 2 5 ij tr(T 2 v i ®v j ) 

+aV{ct) 1 )} + ^[A^dutriTW ® v J ) + aV^ 1 )}}^ 
=: d + eii + fix 2 (2.21) 
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Hence 

det(g)[C - YM C] = af7r 4 + {ae + fb)7f 3 + {ad+fc + eb)7c 2 + {bd + ec)7c + cd . 

(2.22) 

Finally we have thus, using M ij := g 2 B 2 {2K) A /AB i a B{d IJ (E ( }E b j + BfB b j) as 
an abbreviation 

det(g)C 

= {(tr(S 2 - R 2 )) 2 tr(M) - tr(A )tr(M) + 2tr(A M) - 2tr(S 2 - R 2 ) 

tr((S 2 — R 2 — [S, R])M) + cd} 
+ {-Atr(S 2 - R 2 )tr(RT)tr(M) - tr{A x )tr(M) + 2tr(A 1 M) + 2(tr(S 2 - R 2 ) 

tr((RT + TR- [S, T])M) + 2tr(RT)tr((S 2 - R 2 - [S, R])M)) + bd + ec}ir 
+ {-2(tr(S 2 - R 2 )tr(T 2 ) - 2(tr(RT)) 2 )tr(M) - tr(A 2 )tr(M) + 2tr(A 2 M) 

-2(tr(S 2 - R 2 )tr(T 2 M) + tr(T 2 )tr((S 2 - R 2 - [S, R])M) 

-2tr{RT)tr{{S 2 - R 2 - [S, R))M)) +ad + fc+ eb}n 2 
+ {Atr(RT)tr(T 2 )tr(M) - tr(A 3 )tr(M) + 2tr(A 3 M) 

-2(2tr(RT)tr(T 2 M) + tr{T 2 )tr({S 2 - R 2 - [S, R])M)) + ae + /6}tt 3 
+ {(tr(T 2 ))Hr(M) - tr(A 4 )tr(M) + 2tr{A A M) - 2tr(T 2 )tr(T 2 M)}n 4 
=: {(tr(T 2 )ftr(M) - tr(A A )tr(M) + 2tr(A 4 M) - 2tr(T 2 )tr(T 2 M)} 

[vr 4 + avr 3 + (3ti 2 + 7 7r + 5] . (2.23) 

The coefficients a, (3, 7, 5 can now directly be plugged into the resolution for- 
mulas for a general quartic equation due to Cardano and Ferrari as outlined 
in the appendix. 

We have thus arrived at the general solution of the initial value constraints for 
arbitrary matter coupling and nondegenerate 3-metrics and magnetic fields 
as far as the vector and scalar constraint are concerned. The matrix A and 
the field 7r are expressed in terms of the other fields. Note that one cannot 
solve the scalar constraint for one of the 6 independent components of the 
matrix S which appears in fifth order due to the term cd in (2.22) (according 
to Galois theory (see ||) the general algebraic equation of order larger than 
4 is not solvable by radicals). 

In order to complete the solution of the iv constraints one would like to solve 
the Gauss-constraint in terms of the matrix S which however becomes now 
highly nontrivial and is beyond the scope of the present paper. 
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In the appendix we give the above solution for special cases i.e. when the 
Yang-Mills field is absent. In particular, if one only couples the Klein-Gordon 
field, then the scalar constraint depends only on it 2 and the above formulas 
simplify tremendously. 

3 Solving all constraints 

3.1 Solution by using the scalar fields 

First of all we show how a Higgs-field can be used for any semi-simple gauge 
group with rank i to satisfy the Yang-Mills gauge constraint (note that for 
semisimple groups /jjk := fu L diK is totally skew; see ref. [0] for the 
necessary Lie algebra terminology). 
The YM Gauss constraint is 

g I :=V a E^ + f IJ K (f ) J 7i K = 0. (3.1) 

Let V be the complex vector space of dimension dim(G) which is the rep- 
resentation space of the adjoint representation of LG. Let further Y : = 
J T/, (Tj) JK = -f IJK and V ± ((j) 1 ) = {X G LG; ad(Y)X := [Y,X] = 0} 
which shows that V L depends only on <fr (even the dimension k of V 1 - de- 
pends, in general, on the specific element Y : choose a Cartan subalgebra 7i 
of LG and choose the corresponding Weyl canonical form. If Y G 7i then 
1 < dim(V ) = £, if Y is a nonzero-root vector relative to 7i then dimiV^) 
is not characterizable by i only but will depend on the algebra and the root 
chosen. One only knows then that dim(V^) > £ + 1, see ref. 0, where 
VW :=V -V L ). ' 

We will regard the fields as Lie algebra valued by the identification $ l Ti =: <fi 
etc. 

Contracting equation (3.1) with X 1 G V yields 

X T V a Ef = (3.2) 

because X 1 fij<fr kk = tr(Tx[X, (j)])7T K = tr{ir\X, <fi]) = for any it. Hence 
the internal- vector density T> a Ej is weakly 'orthogonal' to the subspace V . 
Hence, the part of the Gauss-constraint which is 'parallel' to V ± has to be 
satisfied independently of the momentum 717. We will satisfy the constraint 
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eqns. (3.2) which can be read as a condition on Ef by employing the following 

new trick which is at the heart of the present approach : 

Let 

Ef := e abc V b v cI (3.3) 

where the generalized tensor v a i is yet arbitrary. Then, using the torsion- 
freeness of the (purely metric-determined) Riemann-connection T a be which 
acts on tensor indices only (the torsion due to the spinor fields shows up in uj l a , 
more precisely in the spin-connection, and is part of the reality conditions, 
see refs. 0,0) we have 

X T V a E^ = g 2 f u K X I B a3 v aK = VX G V ± (3.4) 

which is a purely algebraic restriction for v a i and which can be satisfied iden- 
tically by solving for k components of v a j in terms of the others and of the 
magnetic as well as the Higgs fields. We assume that this has been done in 
the following. The electric fields remain, however, independent of the Higgs- 
momenta which is important for the sequel. 

Are there any restricions implied on the Lie- algebra valued 2-form £ a &j := 
Ej€ ca b by representing it as a exterior differental of the Lie-algebra valued 
1-form v a i i.e. tj = T> A v\ ? The purely algebraic properties of both gener- 
alized tensors are the same, however one has to worry about the 'generalized 
integrability conditions' obtained by taking the exterior differential of the 
last equation 

V Ati — -fj JK Fj A v K . 

The latter equation can now only be satisfied for an arbitrary ti if the mag- 
netic fields Fj are non-degenerate. This restriction does not directly show up 
in equation (3.3) because even for Aj = the rhs of eq. (3.3) is nonvanishing 
for a suitable choice of vj. Hence we obtain the same restriction as CD J on 
the magnetic fields in order that the new trick works, at least when there is 
a non-trivial Yang-Mills-potential. 

We decompose the dual internal vector density 717 into a part 'parallel' and 
'orthogonal' to V 1 - 

7Tj := nj- + vr| , (3.5) 

where ti 1 - G V ± , 7r" G are dual internal vector densities. Hence it follows 
when inserting into the Gauss-constraint 

V a Ef + Y/4 = , (3.6) 
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where it is understood that Ef is replaced by equation (3.3). Equation (3.6) 
can be solved for 7r" since Y is regular on V" while ir 1 - remains unspecified. 
Hence we can regard the gravitational and YM-Gauss law as identically sat- 
isfied in terms of 2 of the gravitational and dim(G)-k components of the YM 
Higgs field and 1 component of the gravitational and k components of the 
YM electric field. Note, however, that there is the additional spin-density 
7r T Tiifj contained in the gravitational Gauss constraint which upon solving 
the Gauss constraint becomes part of the vector v a %- 

When inserting these solutions into the vector constraint, one obtains again 
a genuinely differential relation between the relevant momenta v a i, v a i and 
thus the vector constraint fails to be algebraically solvable in terms of the 
gravitational field. However, we can make use of the gravitational field 
and/or the YM field irj- and 3 Klein-Gordon momenta ir a a — 1, 2, 3 in order 
to solve the rest of the constraints by algebraic methods (the usage of scalar 
fields in order to solve all constraints of GR plus scalar matter is most con- 
venient in the old ADM variables). One solves first the 3 vector constraints 
in terms of 7T a which will then depend linearly on the fields tt^-, Hence, 
the scalar constraint also depends only quadratically on ir^-, 7ij~. Accord- 
ingly, by coupling suitable matter, it turns out to be possible to solve all the 
constraints explicitely. So matter helps to solve the complete iv constraints 
with this method. 

We will now give the explicit formulas. Let ir^ =: h E (f)i and |0| 2 := 
The vector constraint reads 

V a = [ni b P? + F^E*} + njV^ 1 + -Le^iVtP^V^) 

+ir T V a ^} + hE^Vaft + 7r a <j) a a (3.7) 

where it is understood that (f>r,Ef, P® are expressed in terms of the other 
fields as derived above, the crucial point being that they do not depend on 
h,E,TTa- We write this as a matrix relation 

u + h E u + P(7ii,n 2 ,7i 3 ) T = (3.8) 

where the matrix P, which consists of the covariant derivatives of the fields 
<p a , is in general non-singular so that equation (3.8) can be inverted to give 

7T Q = -(P-y a [u Q + hu\a =: (a + bh E ) a . (3.9) 
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Here the coefficients a a , b a do not depend on h E . 

We finally insert this into the scalar constraint (again it is to be understood 
that one has to insert the above solutions for P", Ef, (J) 1 and that q ab is ex- 
pressed in terms of P") and obtain 

C = {-^P/P^-^ 

-A{kP^ t t{D^ + det(q)V^ A , n A )) + ^-q ab d IJ (EjEj + B?B b j) 

+ll(^ijk<f> j (V b P b k )) 2 + det(q)(5 ljq ab (V a ^)(V b ^) + Vtf))] 

+l[d IJ nmj + &et{q){d LJ q ab {V a ^){V b ^ J ) + V^ 1 ))]} 

+ {a a ba}hE + {^b 2 a + ^\(P\ 2 }(h E ) 2 

= : a(h E ) 2 + f3h E + 1 (3.10) 

which yields a quadratic equation for h E - Hence, the inclusion of a gravi- 
tational and YM Higgs-field together with this new method enables one to 
obtain the general solution of the initial value constraints for arbitrary matter 
coupling. Note that the new method introduced is completely independent 
of the CDJ-framework since the CDJ-matrix does not enter the game at any 
stage. 



3.2 Solution without using the scalar fields 

The method of section 3.1 (as well as of section 2) has the unattractive feature 
that explicit use was made of the (possibly spurious) scalar fields. We now 
show that one can do without them altogether by essentially the same trick 
if one uses Weyl fields to solve the vector and scalar constraint as well as the 
gravitational Gauss constraint. 

Inserting the basic ansatz (3.3) into the Gauss constraint yields the purely 
algebraic relation 

V a E1 = g 2 f u K B aJ v aK = -fjj k <P j it k VX e V 1 (3.11) 

for v a i- There are 3 x dim(G) independent components of the electric field 
contained in v a i and dim(G) constraints, so we can choose dim(G) of the 
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components of v a j to depend on the others as well as on the magnetic fields, 
the Higgs field and the V" part of the Higgs-momenta. We assume that this 
has been done from now on. 

Now one inserts formula (3.11) into the rest of the constraints. We assume 
that at least 4 spinor fields are present and solve these remaining 7 con- 
straints in terms of 7 components of spinor field momenta by using similar 
algebraic methods as in eqn. (3.8) ff. If the spinor momenta enter the scalar 
constraint only linearly, i.e. V(v/j a , tta) = V(^ A ), there is not even the need 
to solve a quadratic equation which simplifies the relevant formulas tremen- 
dously. Hence, one could then write the vector and the scalar constraint as 
one matrix relation as in eqn. (3.1.8) for the 7 components of spinor mo- 
menta chosen and solve for them by methods of linear algebra. If this is not 
the case, then we require that the spinor potential is of at most 4th order 
in the spinor momenta (without spatial derivatives) in order that the scalar 
constraint be solvable by algebraic methods. 

We refrain from giving the explicit formulas because the procedure how to 
get them is identically the same as the one of section (3.1). 

Note : in case of the gravitational Gauss constraint one can interprete the 
trick, eqn. (3.3), in the following nice way : 

Provided that the magnetic fields are non-degenerate we can write v a % = 
m a bB^, m a b being an arbitrary tensor density of weight -1. Then the Gauss- 
law yields 

V a P t a = e tjk B]B b k m ah = BBy hc m ab = -e^0V fe (3.12) 

such that the antisymmetric part of the tensor density m is expressible in 
terms of the gravitational Higgs-field 

e abc m ab = -ii^ fe 0V fe (3.13) 

while its symmetric part remains unspecified. Unfortunately, the symmetric 
part enters the vector and scalar constraint differentiated which excludes the 
possibility to solve for m^, at least by algebraic methods. 
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A Appendix 

A.l The resolution formula for quart ic equations 

For the benefit of the reader we include here a brief recipe for solving the 
general quartic equation. The formulas were first found by Cardano and 
Ferrari (see 0). 

1) Turn the general form of a quartic equation 

x 4 + ax 3 + f3x 2 + >yx + 5 = (A.l) 
into its normal form by substituting x = y — a/4 : 

y 4 + py 2 + qy + r = (A. 2) 

2) Make the ansatz 

y 4 + py 2 + qy + r = (y 2 + P) 2 - {Qy + R) 2 (A.3) 

This factorizes into a product of two quadratic equations and can be solved 
by standard methods. Comparison of coefficients yields 

p = 2P -Q 2 , q= -2QR, r = P 2 -R 2 (A.4) 

and results in the so-called cubic resolvent 

p3 _ p2 _ r p + q 2 + _ rp = . p3 + a p2 + bp + c=0 ^ 

2 8 2 

3) Solve the cubic equation which upon substitution P = t — a/3 adopts its 
normal form t 3 + et + / = 0. Let D := (e/3) 3 + (//2) 2 (the discriminant) 
and 

u := f2 + VD, v := ^-{~^D (A.6) 
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then the 3 roots are given by 

U + V . r—U — V U + V . r-U — V - . * 

ti = u + v,t 2 = h tV3 - , t 3 = - iV3 - (A. 7) 

4) Use one of the roots t 1 ,t 2 ,t 3 of the cubic resolvent to determine P,Q,R 
and proceed with formula (A. 3). 

Since already a, j3, 7, <5 look horrible when expressed in terms of the inde- 
pendent components of the fields (compare formula (2.23)) we refrain from 
giving x in terms of a, j3, 7, 5 explicitely and rather discuss a feasible example. 

A. 2 The extended method of CD J applied to a special 
case 

First we switch off only the Yang-Mills field (and, necessarily, its associated 
Higgs-field) . Then the scalar constraint reduces to (recall formula (2.21)) 

fir 2 + 67T + d = (A.8) 

where Ej = A ! a = n 1 = (p 1 = 0, also in the expressions for a,b,c - compare 
formula (2.18). This is now only a quadratic equation for n. 
Although the degree of the scalar constraint cannot be lowered further with- 
out switching off the Klein-Gordon field altogether, a tremendous simplifica- 
tion occurs when one retains only the gravitational and the KG field because 
then the matrix R vanishes. We obtain 

det(^) = [det(S)] + n 2 [tr(ST 2 ) - Ur{S)tr{T 2 )] =: (arc 2 + c)/{2ikB) (A.9) 

i.e. b=0. The scalar constraint reduces to 

C = -KB((tr(S)) 2 - tr(S 2 ) - n 2 tr(T 2 )) 

+ i[vr 2 - AK 2 tr(S 2 v ®v) + (an 2 + c)V((j>))] 

= {-KB(((tr(S)) 2 - tr(S 2 )) + -\-AK 2 tr(S 2 v ® v) + cV (<!>)]} 

2 

+{KBtr(T 2 ) + -[1 + aV((f))]}n 2 

= fn 2 + d (A.10) 

i.e. e=0, the momentum ir enters the scalar constraint without a linear term 
! The matrix T is simply = B%/ ' {2B)(f)^ a eijk- 
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